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Abstract 

In the Next-to-Minimal Super symmetric Standard Model, neutralino dark matter 
can annihilate into a pair of photons through the exchange of a CP-odd Higgs boson 
in the s-channel. The CP-odd Higgs boson couples to two photons through a loop of 
dominantly higgsino-like charginos. We show that the parameter space of the NMSSM 
can accommodate simultaneously i) neutralino-like dark matter of a mass of about 
130 GeV giving rise to a 130 GeV photon line; ii) an annihilation cross section of or 
larger than 10 -27 cm 3 s _1 ; iii) a relic density in agreement with WMAP constraints; 

iv) a direct detection cross section compatible with bounds from XENON100, and 

v) a Standard Model like Higgs mass of about 125 GeV. However, the CP-odd Higgs 
mass has to lie accidentally close to 260 GeV. 



1 Introduction 



The fact that annihilation of dark matter can give rise to photons with a well-defined 
energy [1-22] has motivated searches for such gamma lines originating from the center of 
our galaxy by the EGRET [23], H.E.S.S. [24,25] and Fermi LAT [26-28] experiments. 

Recent analyses of the publicly available data from Fermi LAT [29-32] have discovered 
hints for a gamma line at E 1 ~ 130 GeV in the form of an excess of about 3 — 4 standard 
deviations, assuming that the background flux can be approximated by a single power law. 
An interpretation of this excess as dark matter pair annihilation into a pair of photons 
would require a partial annihilation cross section of about 10 -27 cm 3 s _1 [30,31]. 

However, more general parametrizations of the background flux [33, 34] reduce the sig- 
nificance of the excess, making it compatible with a diffuse background possibly of instru- 
mental or astrophysical origin. The Fermi LAT collaboration preferred to interpret the 
same data only in terms of upper bounds on the partial annihilation cross section [28]. 

Still, future additional data could confirm the present hints for a possible excess; hence 
it is of interest to study whether it could be explained within concrete models for dark 
matter which are compatible with bounds on its relic density from WMAP [35] and bounds 
on its direct detection cross section (in the relevant mass range) from XENON 100 [36]. 

Following the publication [30], different types of such models have been proposed: mod- 
els with an extra U(l)' gauge symmetry where the dark matter couples only to the extra 
Z' gauge boson, and the Z' to photons via a Chern-Simons term [37]; models with an extra 
singlet and extra charged fields allowing for an enhanced branching ratio of the Standard 
Model-like Higgs boson into two photons [38]; extensions of the Minimal Supersymmet- 
ric Standard Model (MSSM) by right-handed neutrinos and extra Higgs doublets (with a 
right-handed sneutrino as dark matter) [39]; decaying dark matter in the MSSM extended 
by additional fields and couplings [40]; extensions of the Standard Model or the MSSM by 
singlets with Peccei-Quinn symmetry where dark matter can annihilate into a photon pair 
via an axion in the s-channel [41]; a string/M-theory motivated version of the MSSM (with 
an unconventional spectrum involving very heavy scalar superpartners) where wino-like 
dark matter can annihilate into a photon plus a Z boson [42]. (In the MSSM with bino- or 
higgsino-like dark matter, the annihilation cross section into one or two photons would be 
too small [2-5, 7].) 

In the present paper we consider the simplest supersymmetric extension of the Stan- 
dard Model with a scale invariant superpotential, the Next-to-Minimal Supersymmetric 
Standard Model (NMSSM) [43,44]. In the NMSSM, the Higgs/higgsino mass term piH u H d 
in the superpotential of the MSSM is replaced by a coupling \SH u Hd to a gauge singlet 
superfield S. Compared to the MSSM, the Higgs sector of the NMSSM includes additional 
neutral CP-even and CP-odd states. It is known that a large dark matter (neutralino) 
annihilation cross section into two photons can arise from neutralino annihilation through 
the NMSSM-specific CP-odd Higgs boson in the s-channel [9,20,21,41]. Moreover, due 
to the coupling XSH u Hd the lightest CP-even Higgs mass is naturally heavier than in the 
MSSM [43-48] which makes it easier to explain a Higgs mass of 125 GeV [49-59]. 

Hence an obvious question is whether the parameter space of the NMSSM can describe 
simultaneously the following phenomena: 1) A neutralino-like dark matter particle of a 
mass of about 130 GeV with a partial annihilation cross section into two photons of about 
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10 -27 cm 3 s _1 and a relic density compatible with WMAP constraints, 2) a Standard 
Model-like Higgs boson with a mass of about 125 GeV, but 3) respecting the upper bound 
of XENON100 on the dark matter direct detection cross section. 

We find that this is indeed possible, provided the NMSSM-specific coupling A is relatively 
large, A ~ 0.6. This coupling plays several roles simultaneously: i) it determines the 
coupling of the neutralino-like dark matter to the NMSSM-specific CP-odd Higgs boson; ii) 
it determines the coupling of the CP-odd Higgs boson to charged higgsinos, whose loop leads 
to a large decay width into two photons; iii) it helps to increase the mass of the Standard 
Model-like CP-even Higgs boson. However, the mass of the NMSSM-specific CP-odd Higgs 
boson has to be close to 260 GeV (the s-channel pole) for a neutralino annihilation cross 
section into two photons equal to or larger than 10~ 27 cm 3 s -1 (see [62]). On the other 
hand, no new particles or interactions beyond the NMSSM need to be introduced if the 
present hints for a 130 GeV gamma line from dark matter annihilation in our galaxy are 
confirmed. 

In the next section we introduce the NMSSM and describe, which properties of the 
NMSSM allow to describe simultaneously all the above phenomena. We indicate a possible 
region in the parameter space, and discuss the range of possible direct detection cross 
sections and neutralino annihilation cross sections into two photons. Details of the Higgs 
and neutralino sector are presented for a typical benchmark point. Section 3 is devoted to 
a short summary and conclusions. 

2 A 130 GeV Photon Line in the NMSSM 

The NMSSM differs from the MSSM due to the presence of the gauge singlet superfield 
S. In the simplest realisation of the NMSSM, the Higgs mass term [iH u H d in the MSSM 
superpotential Wmssm is replaced by the coupling A of S to H u and H d , and a self-coupling 
kS 3 . Hence, in this version the superpotential H^nmssm is scale invariant, and given by: 

Hnmssm = ^SH U ■ H d + — S 3 + . . . , (1) 

where the dots denote the Yukawa couplings of H u and H d to the quarks and leptons as in 
the MSSM. The NMSSM-specific soft SUSY breaking terms consist of a mass term for the 
scalar components of S, and trilinear interactions associated to the terms in Wnmssm : 

-4mssm = ™II<S1 2 + (\A x H u -H d S+^A K S 3 ) + h.c. . (2) 
Subsequently we define the vacuum expectation values (vevs) 

(H u ) = v u , (H d )=v d , (S) = s. (3) 
In terms of s, the first term in W^nmssm generates an effective //-term with 

Aieff = As . (4) 
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Using the minimization equations of the potential in order eliminate the soft SUSY 
breaking Higgs mass terms, the Higgs sector of the NMSSM is characterized (at tree level) 
by the six parameters 

A , k , A x , A K , /i eff , tan (3 = — . (5) 

v d 

The neutral CP-even Higgs sector contains 3 states Hi, which are mixtures of the CP- 
even components of the superfields H u , H d and S. Their masses are described by a 3 x 3 
mass matrix Ai 2 Hi j, where the dominant contribution to the singlet-like component -M^ 33 
reads [43,44] 

M H33 ~ ks(A k + 4ks) . (6) 



The neutral CP-odd Higgs sector contains 2 physical states A t , whose masses are de- 
scribed by a 2 x 2 mass matrix M. Ai j where M.\u corresponds to the MSSM-like CP-odd 
Higgs mass squared. The dominant contributions to the singlet-like component M. 2 A22 an d 
the singlet-doublet mixing term M. 2 A12 are given by 



M 2 A22 ~ -3ksA k , M 2 Al2 ~ X(A X - 2ks)^vI + v 2 , (7) 
respectively [43,44]. 

In the neutralino sector we have 5 states Xii which are mixtures of the bino B, the 
neutral wino VU 3 , the neutral higgsinos H d , H u from the superfields and H®, and 
the singlino from the superfield S. Their masses are described by a symmetric 5x5 mass 
matrix -M^o^ given by 



M x 0ij 



( M 1 -af aj^ ^ 

M 2 -eg- 

-yU cff -\V U 

-Xv d 

V 2ks ) 



(8) 



where Mi, M 2 are the soft SUSY breaking bino and wino mass terms. The lightest eigenstate 
Xi of the neutralino mass matrix (8) is considered as the dark matter particle. In the above 
basis, its decomposition is written as 



X °i = N n B + Ni 2 W 3 + Ni 3 H d ° + N U H U ° + Ni 5 S . (9) 

The pair annihilation of Xi into two photons can be generated via a dominantly singlet- 
like CP-odd Higgs state As in the s-channel [9,20,21], see Fig. 1. (Similar diagrams with 
X ± replaced by (s)quarks or (s)leptons are numerically negligible.) 

Now we turn to the masses and couplings of the particles and the corresponding regions 
in the NMSSM parameter space, which have all desired phenomenological properties. First 
we consider the dominantly singlet-like CP-odd Higgs state A$. For a sufficiently large Xi 
pair annihilation cross section into two photons (and = M x o ~ 130 GeV), we need 

M As ~ 2M X? ~ 260 GeV , (10) 



3 



Figure 1: The dominant diagram for pair annihilation of xl into two photons via a mostly 
singlet-like CP-odd Higgs A s in the NMSSM. 

which can be achieved by appropriate values of — 3ksA k . Moreover, the SU(2) doublet 
admixture of As must be small: Otherwise tree level diagrams similar to Fig. 1 but with 
As decaying directly into bb (or into Z plus a light CP-even Higgs boson), lead to a too 
large pair annihilation cross section of Xi such that its relic density is below the WMAP 
bound. From the second of Eqs. (7) the mixing of As with the MSSM-like doublet is small 
for 

A x ^2ks. (11) 

Next we consider the dark matter particle Xv It would have a large singlino component 
for small 2ns. However, from Eq. (6) and the first of Eqs. (7) one can derive 

(2ks) 2 ~ M 2 H33 + l -M\ s ■ (12) 

from Ma s ~ 260 GeV and Ai 2 H33 > it follows that 2ks cannot be small. Hence, assuming 
Mi < M 2 /2 (assuming universal gaugino masses at the GUT scale), it follows that Xi 
has dominant bino and higgsino components. A priori large higgsino components seem 
desirable, given the required coupling of Xi to As in Fig. 1: This coupling is induced 
by the first term XSH U ■ H d in the superpotential (1), which leads to a Yukawa coupling 
\AsH u Hd . Likewise, the coupling of As to the charginos x ± originates from the higgsino 
components H u + , H d of x ± an d the Yukawa coupling \A S H U + H d . 

However, too large higgsino components of Xi imply again a too small relic density; 
diagrams with charginos and neutralinos in the t-channel (and W + W~ or ZZ in the final 
state), CP-even Higgs bosons in the s-channel etc. would lead to a too large pair annihilation 
cross section of Xv Hence we end up with a dominantly bino-like Xi> but with non-zero 
(non-negligible) higgsino components. Its mass of 130 GeV has to follow from appropriate 
values of Mi and fi c s, with M 1 < /^ e ff. 

Finally we require a SM-like CP-even Higgs boson Hsm with a mass Mh sm near 
125 GeV. Although its existence is not confirmed at present, it is interesting to inves- 
tigate whether it could comply with the above properties of the neutralino and CP-odd 
Higgs sector. It has been known since a long time that the SM-like CP-even Higgs boson 
can be heavier in the NMSSM compared to the MSSM due to the NMSSM-specific coupling 
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\SH u H d [43-48], provided A is large and tan/3 is relatively small. Subsequently we choose 
corresponding values of A and tan/3 such that M HsM ~ 125 GeV [49-59]. We find that the 
above properties in the neutralino and CP-odd Higgs sector imply that H S m is the lightest 
CP-even Higgs state; the singlet-like CP-even Higgs state has a mass > 200 GeV. The 
scenario with a lightest singlet-like Higgs state and a next-to-lightest SM-like Higgs state 
(allowing for an enhanced branching ratio into 77 [51,53,55,56]) seems difficult to realize. 
Herewith we have sketched the interesting regions in the parameter space in (5). 

An open question remains whether the dark matter particle can comply with the con- 
straints from XENON100 on its spin-independent detection cross section: x^-nucleon scat- 
tering is induced dominantly by if^M-exchange in the t-channel [67], and the Hsm~Xi~X\ 
vertex is proportional to the product of the bino- and higgsino-components of Xi (from gi x 
bino x Higgs x higgsino terms in the Lagrangian, where g\ is the U(l)g gauge coupling). 
This issue will be studied below. 

We have scanned the parameter space of the general NMSSM with help of the NMSSM- 
Tools package [68,69], supplemented with suitably modified formulas for the cross sections 
f° r x\x\ ~~ >" 77 (and for x\x\ ~^ %l) from [2,4]. MicrOmegas [70-72] is used for the cal- 
culation of the dark matter relic density and direct detection cross sections. For the latter 
we have to specify the strange quark content of the nucleons, i.e. the relevant sigma terms. 
We use the most recent values from [73] with, to be conservative, a value for a^N near the 
lower end of the la error bars: a„N = 40 MeV, o"o = 39 MeV. 

For the soft SUSY breaking terms we made the following choices: 

• Squark masses of 1.5 TeV, except for the left-handed squarks of the 3rd generation 
(1 TeV) and the right-handed top squark (300 GeV). The latter values are motivated 
by universal soft scalar masses at the GUT scale [57], and alleviate LHC constraints 
from direct SUSY searches due to the more complicated squark and gluino decay 
cascades involving the light stops. 

• Trilinear soft susy breaking terms At = Aj, = — 1.1 TeV. 

• Slepton masses in the 140 — 500 GeV range such that the SUSY contributions to 
the anomalous muon magnetic moment are sufficiently large (inspite of low values of 
tan/3), while slepton exchange in the t-channel of the pair annihilation cross section 
of Xi does not imply a too small relic density. 

• Whereas we vary M\ in the 140 — 160 GeV range (see below), M2 and the gluino mass 
M 3 are kept fixed at M 2 = 300 GeV, M 3 = 800 GeV for simplicity. 

• Finally we use 173.1 GeV for the top quark pole mass. 

Subsequently we impose the following phenomenological constraints: 

• M x o = 129 - 131 GeV and (av) {xiXi ~> 77) > 10~ 27 cm 3 s^ 1 in order to obtain a 
photon line in agreement with the excess found in [30,31]. 
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• Upper bounds on annihilation cross sections into W + W~ , ZZ, bb and rf channels 
from the Fermi LAT collaboration [28,60] (see also [61,62]), as well as bounds from 
PAMELA on the antiproton flux [63]. For the determination of these cross sec- 
tions/fluxes we use micrOMEGAs2.4 [64]. 

• A relic density complying with the WMAP bound Vth 2 = 0.1120 ± 0.011 [35] (with 
2ct error bars). 

• A SM-like Higgs boson with Mh sm = 124 — 127 GeV, as it was confirmed recently by 
the ATLAS and CMS collaborations [65,66]. 

• Constraints from B-physics as implemented in NMSSMTools (which have actually no 
impact for the low values of tan/3 considered here). 

• A sufficiently large SUSY contribution Aa M to the muon anomalous magnetic moment 
as implemented in NMSSMTools. 

• Constraints from the absence of Landau singularities of the running Yukawa couplings 
below the GUT scale, and the absence of unphysical global minima of the Higgs 
potential. 

We have found corresponding points in the NMSSM parameter space, both below and 
slightly above the present XENON100 bound on the spin- independent dark matter - proton 
cross section of a{p)si ^, 1-2 x 10~ 8 pb for M x o ~ 130 GeV [36]. In Fig. 2 we show a(p)si 
as function of Mi for a sample of such points, where we varied the parameters in (5) in the 
range A = 0.6 - 0.615, k = 0.326 - 0.329, A x = 240 - 400 GeV, A K = -130 - (-60) GeV, 
// eff = 230 - 445 GeV, tan/3 = 1.68 - 1.82. 

In Table 1 we show the details (parameters, masses and relevant observables) for a sam- 
ple point with Mi = 150 GeV. The couplings of Hi = H SM to quarks, leptons, electroweak 
gauge bosons, gluons and photons are SM-like within ~ 5%. denotes the coupling of 
As to 6-quarks normalized to the one of a SM-like Higgs boson; its small value underlines 
its singlet-like nature. 

The following remarks are in order: The larger Mi, the smaller one has to choose n e s in 
order to maintain M x o ~ 130 GeV. It follows that, for larger Mi, the higgsino component 
of Xi increases leading to a larger XiXi annihilation cross section. Hence, for too large Mi, 
the relic density falls below the WMAP bound. 

On the other hand, for smaller Mi one has to choose larger values for /x e ff implying a 
smaller higgsino component of Xi, which explains the decrease of cr(jp)si in Fig- 2. However, 
simultaneously the coupling of A s to xi decreases as well. As a consequence, the mass Ma s 
of As has to be closer and closer to the pole 2M x o in order to obtain cr(xiXi ~~ >* 77) > 
10 -27 cm 3 s" 1 . In order to clarify the required tuning, we show (t(XiXi ~^ 77) i R Fig- 3 as 
function of Ma s for the point listed in Table 1. (Note that the finite width T(Ag) ~ 1.6 KeV, 
dominated by As — > 77, is not visible in Fig. 3. Due to the small couplings of As to quarks 
and leptons, the contributions of As in the s-channel to the annihilation cross sections into 
bb and rr final states are well below the Fermi LAT bounds. The annihilation cross sections 
into W + W~ and ZZ originate from the second CP-even Higgs boson in the s-channel. The 
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antiproton flux has a maximum of ~ 2.47 x 10~ 4 (GeVm 2 ssr) -1 for an energy of ~ 2 GeV, 
which is well below the PAMELA bound.) 

We see that &(xiXi ~~ >" ll) is larger than 10~ 27 cm 3 s _1 only within a ~ 0.7 GeV 
wide window of Ma s - This required tuning becomes worse for lower values of Mi, and 
would be the price to pay for a possibly stronger constraint on cr(p)si from XENON100 
in the future. (On the other hand, modifications of the present best estimates for the 
Higgs-nucleon coupling and/or the local dark matter density could alleviate the present 
constraints from XENON100.) 

Finally we should add that diagrams similar to Fig. 1, but with one photon replaced 
by a Z boson, contribute to <t(XiXi ~~ >* %l) leading to an additional photon line with, 
for M x o ~ 130 GeV, ~ 114 GeV. For the present scenario we find <j(xiXi ~^ %l) ~ 
75% x a(xiXi ~^ ll)- Such an additional line would be compatible with the structure 
observed in [32]. In any case, additional lines — also from cr(XiXi ~~ ^ Hj) or interpreting 
the 130 GeV line as due to XiXi ~^ %1 — can be interesting checks of such scenarios in the 
future [74]. 
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Figure 2: u{p)si as function of Mi for a sample of points satisfying all other phenomeno- 
logical constraints as indicated in the text. The horizontal line indicates the present bound 
from XENON100 for a 130 GeV dark matter particle, which holds for all points shown. 
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Parameters 


A 


0.61 


K 


0.328 


A x 


267 


A K 


-114.1 


tan (3 


1.8 




269 


Mi 


150 


left-h. slepton masses 


150 


right-h. slepton masses 


160 


A e = A^ = A T 


500 


Sparticle masses 


ITlg 


971 


( m <j) 


1530 


mi! 


204 


m h 


1034 


m l, 


1005 


m H 


154 


Xi 


129.6 


M v o 

A2 


217 


M x5 


287 


M x2 


309 


M xS 


376 


M xt 


210 


M xf 


370 



Higgs masses 


Mm (= Mm ) 


124.3 


M H2 


256 


M Hz 


519 


M Al (= M As ) 


258.9 




3 x 10~ 3 


M A2 


515 


M H ± 


511 


Components of Xi 


N 2 


0.826 


N 2 


0.026 


N 2 

ly 13 


0.077 


N 2 


0.065 


N 2 


0.009 


Observables 


Qh 2 


0.11 


a(p) SI [10- 8 pb] 


1.21 


(<™)(XiXi -> 77) [10~ 27 cm 3 s" 1 ] 


1.1 


(™)(XiXi -> Zj) [10- 27 cm 3 s- 1 ] 


0.8 


M(XiXi ~* WW ) [10" 27 cm 3 s" 1 ] 


3.46 


M(XiXi ^Z) [10- 27 cm 3 s" 1 ] 


0.26 


M(XiXi -> 66) [10- 27 cm 3 s" 1 ] 


0.60 


(<™)(XiXi -» rf) [10~ 27 cm 3 s" 1 ' 


0.09 


Aa M [10- 10 ] 


6.5 ±3.0 



Table 1: Properties of a sample point with Mi = 150 GeV. Dimensionful parameters are 
given in GeV. R b ^ s denotes the coupling of As to 5-quarks normalized to the one of a SM- 
like Higgs boson. The components of Xi are defined in Eq. (9). The value of Aa M includes 
theoretical error bars. 



3 Conclusions 

In this paper we have shown that the simplest version of the NMSSM (with a scale invariant 
superpotential) could explain a 130 GeV photon line from dark matter annihilation with 
a(XiXi ~^ 77) > 10~ 27 cm 3 s _1 and, simultaneously, a 125 GeV SM-like Higgs boson. No 
additional fields or couplings need to be introduced. All constraints from WMAP on the 
relic density, from XENON100 on the direct detection cross section, from colliders and from 
precision observables can be satisfied. 

However, the mass Ma s of the singlet-like CP-odd Higgs scalar As has to satisfy ac- 
cidentally Ma s ~ 2M x o ~ 260 GeV to a precision < 1 GeV. This "fine-tuning" would 
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Figure 3: <j{xiXi ~* 77) as function of M^ s for the point listed in Table 1. 

become worse if bounds on the direct detection cross section become stronger (but could 
be relaxed otherwise). 

Unfortunately a direct verification of this scenario at colliders through searches for a 
130 GeV photon line seems hopeless: Due to the singlet-like nature of As, production cross 
sections for this state (as well as decay widths of sparticles or other Higgs bosons into this 
state) are too small. Only the mass of 130 GeV of the LSP Xi should fit the data, once 
searches for supersymmetry turn out to be successful. Of course, first of all the present 
hints for a 130 GeV photon line [29-32] need to be confirmed. 
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